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We apply the Pseudo-C; formalism to obtain an unbiased, approximate method for efficient 
simultaneous estimation of several cosmological parameters from large, almost full-sky cosmic 
microwave background data sets. 



1 Introduction 



Within the standard model of cosmology there are about 10 parameters which characterise the 
properties of our Universe. It is one of the key goals of future CMB experiments such as MAP 
and Planck to determine these cosmological parameters to high precision. This undertaking faces 
the challenge that realistic CMB data is necessarily incomplete and noisy. The Galaxy obscures 
roughly a third of the sky and because of the smallness of the anisotropy signal, detector noise 
is not negligible in the analysis. This leads to the computational challenge which was expertly 
described at this meeting in the contribution by BorrilH. 

In this talk we apply the pseudo-Q formalisnll to this problem and show that it can be 
used to develop an approximate form of the likelihood which has several useful properties: 
it is Gaussian and hence easy to apply; it does not suffer from the usual disadvantages of 
Gaussian approximations such as obtaining negative estimates of positive definite quantities; 
it is computationally efficient with memory usage of Npix and number of operations scaling as 
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^pix P^^ likelihood evaluation with a very small pre-factor leading to thousands of likelihood 
evaluations per CPU hour. 
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2 The approximation scheme 



Given a true CMB sky T and an experimental setup and observation strategy (encoded in the 
beam pattern B, the survey geometry W and the noise distribution on the sky W]\fT]\f) we can 
represent the observed temperature anisotropy map as 

r(7) = w{-f) [B * r(7) + WNh)TN{7)] (i) 

This temperature field T can be decomposed into spherical harmonics coefficients 

di^= f dnYiUimi). (2) 
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The notation "J^" denotes integration over the fraction of the sky covered by the survey. These 
combine into the observed power spectrum coefficients which we call pseudo-Q, 



Cl = ^^Er^lm\'. (3) 



Inii we derive the exact statistics of the pseudo-Q, under the assumptions of azimuthal 
survey geometry and noise which is uncorrelated from pixel to pixel and whose amplitude varies 
only from latitude to latitude. The results we derived were still a superb approximation for 
strongly non-azimuthal noise patterns. 

We found that in the case of large sky coverage the Pseudo-C; distributions were nearly 
indistinguishable from Gaussian distributions of the same means and variances as long as / > 100. 
The fact that many of the cosmological parameters are sensitive to the power spectrum at 
precisely these small scales led us to propose the following approximation to the likelihood: 



C{Ci) = n exp 
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(4) 



In this approximation, maximum likelihood estimation has reduced to simple fitting, however 
with the correct means and variances. 

To illustrate, we solve the problem of estimating 3 parameters (f^c and Hq) simultane- 
ously from a sky with 12 x 10^ pixels of which 66% are observed. The response of the experiment 
is modelled as a Gaussian beam of FWHM 12 arcminutes. The noise template is inhomogeneous 
and not azimuthally symmetric with rms amplitude of 124fj,K per 3.5 arcminute pixel. To com- 
pare with the naive Gaussian approach and to show that our method is unbiased, we compute 
maximum (approximate) likelihood estimates from 100 realisations of the sky and plot a repre- 
sentation of the empirical distribution of parameter estimates in three dimensions in Figures |l| 
(naive x^) aiid § (our approach). 

Our estimates are unbiased. The distributions of the estimates are clearly centered on the 
true values. 

We stress that our approach avoids the usual difficulties of Gaussian approximations. For 
example, even though we use the Gaussian pdf, which of course does not exclude negative Ci, 
they are assigned an exceedingly small probability. This is because no attempt is made to 
subtract out the noise contribution from the pseudo-C/ — instead it is modelled consistently 
and the (signal x noise) cross term which is present in each realisation is not allowed to dominate. 
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Figure 1: Multi-parameter estimation of Qq, Qt and Ho using naive fitting. The three panels show three views 
from different directions of the empirical distribution of the parameter estimates. Each sphere represents one bin 
of the three-dimensional distribution. The size and shading of a sphere indicates the number of realisations (out 
of 100 total) which led to parameter estimates within its bin. The true parameter values are at the origin of the 
coordinate axes. It is clearly visible that the distribution is shifted with respect to the true distribution by an 
amount which is inconsistent with the width of the distribution. In other words the true values could be ruled 

out at high significance if this estimate was used. 




Figure 2: Multi-parameter estimation of SIq, Qt and Ho using our approximate likelihood Eq. (^). The three 
panels show three views from different directions of the empirical distribution of the parameter estimates. Each 
sphere represents one bin of the three-dimensional distribution. The size and shading of a sphere indicates the 
number of realisations (out of 100 total) which led to parameter estimates within its bin. The true parameter 
values are at the origin of the coordinate axes. It is clearly visible that the distribution is correctly centered on 

the true values. 



